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Abstract 

A planar force-constant model is developed for longitudinal phonons of 

wurtzite GaN and A1N propagating along the [0001] direction. The proposed 

model is then applied to the study of the phonon modes in hexagonal (0001)- 

GaN/AlN supperlattices in longitudinal polarization. The confinement of the 

superlattice phonon modes is discussed. 
PACS numbers: 68.20.Dj, 68.35. Ja 
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I. INTRODUCTION 



Among the III-V semiconductors the nitride-based semiconductors have received both 
scientific and technological attention over the past years owing to their excellent physical 
properties.0 The fascinating mechanical properties, such as hardness, high melting temper- 
ature, high thermal conductivity, and large bulk moduli make them useful for protective 
coatings. In particular, the electronic properties, characterized by large band gaps and rela- 
tively low dielectric constants, indicate tremendous potential applications in optical devices 
working in the blue and ultraviolet (UV) wavelengths. Recently, high-power blue and green 
light emitting diodes (LEDs) have been fabricated by using InGaN-based multiquantum-well 
structures and are commercially available.! Blue and UV nitride semiconductor laser diodes 
(LDs) have also successfully been demonstrated]! 

Consequently there has been great interest in the fundamental properties, the under- 
standing of which can render the improvement of the materials and device quality possible. 
Heterostructures based on GaN and Gai-^ALjN have become important candidates for op- 
toelectronic materials. Great attention has been paid to the nitride heterostructures and 
layered structures. Short-period superlattices (SL's) have been successfully fabricated by 
advanced growth techniques-^ Phonon spectra of these SL's were investigated by infrared 
(IR) spectroscopy! 

It has been known that GaN and A1N are crystallized in the stable wurtzite structure 
and metastably in the cubic structure. Phonons in cubic GaN/AIN systems have drawn 
interest in recent years. Theoretical studies on phonons in cubic GaN/AIN SL's have been 
reported by several groups. There have been, however, few studies of phonons in hexagonal 
GaN/AIN SL's,ll2l where the constituents have wurtzite structures. 

In this paper we develop a planar force-constant model (PFCM) to describe the lon- 
gitudinal phonons in wurtzite GaN and A1N propagating along the [0001] direction. We 
then apply the PFCM to study the phonons in hexagonal (0001)-GaN/AlN SL's. The paper 
is organized as follows. The PFCM is described in Section II. The calculated results for 
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hexagonal (0001)-GaN/AlN SL's are presented in Section III and conclusions are given in 
Section IV. 



II. PLANAR FORCE-CONSTANT MODEL 



In general, a three-dimension calculation is necessary in the discussion of phonons because 
of the coupling between longitudinal and transverse vibrations. For the vibrations of wurtzite 
compounds along the [0001] direction, however, the longitudinal and transverse vibrations 
can be decoupled. As a result, the longitudinal and transverse vibrations can be treated 
separately. 

The PFTM for wurtzite GaN or A1N is schematically displayed in Fig. [|. The dynamical 
matrix elements are given by 



where q is the wave- vector; / is the index of a unit cell and k stands for the index of atoms in 
the unit cell; M K is the mass of the atom k; <f> (Ik; I'k') is the planar force- const ant between 
the atoms (Ik) and (I'k'); U(Ik) is the position of the atom (Ik). The dynamical matrix 
subjects to the following secular equation 



The eigenvalues and eigenvectors of the vibrational modes can be obtained by solving Eq. 



Four planar force-constants shown in Fig. [I] are used for either GaN or A1N. They are 
two cation-anion planar force- const ants, one cation-cation force- const ant and one anion- 
anion force-constant. The on-site force- const ants <P(Ik; Ik) can be obtained from the sum 
rules 
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(2) 



1 




(3) 



3 



The values of the four planar force- const ants are obtained by a fitting procedure. The 
existing experimental data that are useful in fitting are the Ai(LO) phonon frequency at the 
zone-center determined by Raman spectroscopy. Therefore, the four planar force-constants 
are obtained by fitting them to the three-dimensional calculations.0 The fitted Ai(LO) 
phonon frequency, two B x phonon frequencies and the elastic constant C33 are given in 
Table | together with the corresponding experimental data and other theoretical results. It 
is obvious that the fitted Ai(LO) phonon frequencies of both GaN and A1N are in good 
agreement with the Raman experimental data. The two Bi phonon frequencies agree fairly 
well with other theoretical calculation. The fitted elastic constants C33 are also in good 
agreement with the experiment. The resulting four planar force-constants are given in Table 

E 

The calculated phonon dispersion based on the obtained planar force- const ants for 
wurtzite GaN and A1N for the longitudinal polarization along the [0001] direction are shown 
in Fig. 0. It is found that the four planar-force-constants can give very satisfactory phonon 
dispersion of both GaN and A1N compared with other three-dimensional calculations.!^ 



III. CALCULATED RESULTS 

From Fig. |] some general features of phonon modes in hexagonal GaN/AIN SL's can 
be postulated. If a SL mode is located inside GaN or A1N continuum while GaN and A1N 
phonon continua do not overlap, vibrations of this mode will sharply be confined to one of 
the constituent layers since the other constituent cannot sustain this kind of vibrations. If 
a mode is located in a region where GaN and A1N phonon continua overlap the mode will 
be a resonant, quasi- confined one because the vibrations in another constituent layers will 
be excited. If there is a SL mode with its frequency lies neither in GaN continuum nor A1N 
continuum the vibrations of this mode will be constrained in the interface region. It is a 
microscopic interface mode. 

With the PFCM introduced in the above section phonons in any hexagonal (0001)- 
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(GaN) m /(AlN) n SL's can be studied, where m and n are the number of bilayerslB of GaN 
and A1N, respectively. For GaN/AIN SL's, the planar force- const ants are assumed to have 
the same values as corresponding bulk except that across an interface the planar force- 
constants are taken to be the average of those for GaN and A1N. 

The calculated phonon dispersion and vibrational patterns is given in Fig. [3] for a (0001)- 
(GaN)s/(AlN)5 SL. It is obvious that the vibrations of modes 1-3 are sharply confined to 
A1N layers. They are AIN-like modes originated from the Ai branch of A1N. There is nearly 
no dispersion for these modes due to the strong confinement of vibrations. Mode 4 is also 
an AIN-like mode, but it is from the optical Bi branch of A1N. The frequency of mode 
5 is slightly above the Ai branch of GaN. It is an AIN-like mode, from the Bx branch of 
A1N. Finite excitations in GaN layers exist near the interfaces. Modes 6-9 are resonant, 
quasi-confined GaN-like modes due to the fact that their frequencies are inside the overlap 
region of the Bi branch of A1N and Ai branch of GaN. Vibrations in A1N layers could be 
excited. Consequently, there exists finite dispersion for these modes. Mode 10 is a GaN-like 
mode with vibrations confined to GaN layers. Modes 11-13 are AIN-like modes and they 
are from the acoustic Bx branch of A1N. The remaining modes are folded acoustic modes. 
The minigaps in the zone-center and zone-boundary are due to the difference of the elastic 
properties of GaN and A1N. 

No interface modes are found in hexagonal GaN/AIN SL's. This may stem from the fact 
that there are no new bonds across the interface. 

For a confined mode the vibrations are sharply confined to one of constituent layers. 
In consequence, the confined mode can be viewed as the standing wave, described by an 
effective bulk wave-vector 

where n is the nominal number of bilayers, j is the order of a confined mode, do is the 
nearest-neighbor cation-cation or anion-anion distance and S is a parameter describing the 
degree of penetration of the vibrations of a confined mode into the adjacent constituent. 
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The frequency of the confined mode can be derived simply from the bulk dispersion 

^sl = w(qj) (5) 

instead solving the secular equation about the dynamical matrix if the parameter 5 is prop- 
erly chosen. The above discussions are only valid for the confined modes not for the resonant, 
quasi- confined modes or the extended modes. 

To test the above idea we give in Fig. [| the comparison of the bulk dispersion and the 
frequencies of the confined modes of a (GaN) 10 / (A1N) 10 SL mapped according to Eq. (|J). A 
good correspondence is established between the confined modes and the bulk dispersions for 
both GaN and A1N when 5 is chosen to be 1 for all branches. This means that the vibrations 
of GaN or A1N confined modes extend about one bilayer into the adjacent constituent. For 
those resonant, quasi-confined modes significant vibrations exist in the other constituent 
layers. Their frequencies are somewhat dependent on the thickness of the constituents. 
Therefore standing wave picture is no longer valid. This can be clearly seen from Fig. f|. 

IV. CONCLUSIONS 

In the present work a PFCM is developed for wurtzite GaN and A1N for the longitudinal 
phonons propagating along the [0001] direction. We apply the PFCM to the study of the 
phonons in hexagonal (0001)-GaN/AlN SL's. There are confined modes and resonant, quasi- 
confined modes depending on their frequencies. No interface modes are found in this system. 
The vibrational confinement is investigated. For confined modes either GaN-like or A1N- 
like their vibrations can be viewed as standing waves described by an effective wave-vector. 
Their frequencies can be obtained directly from the bulk dispersion. 
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TABLES 

TABLE I. Comparison of the fitted phonon frequencies at the zone-center and elastic constants 
of wurtzite GaN and A1N with other experimental and theoretical results. The phonon frequencies 
are in units of cm -1 and elastic constants are in units of GPa. 







w(Ai) 


"(B?>) 


"(B?>) 


C33 


GaN 


This work 


739 


690 


330 


399 




Others 


735 a , 736 b 


677 c , 697 d 


330 c , 335 d 


398 e 


A1N 


This work 


893 


718 


542 


390 




Others 


893 f , 891 g 


717 c , 703 d 


553 c , 534 d 


389 h 



a Raman measurement, from Ref. 13 



b Raman measurement, from Ref. 114 



Tirst-principles full-potential LMTO calculations, from Ref. [L5 
d First-principles pseudopotential calculations, from Ref. [16. 
e Measurement, from Ref. [T7|. 
f Raman measurement, from Ref. ITS 



g Raman measurement, from Ref. |I9 
h Measurement , from Ref. 



TABLE II. Values of the planar force-constants obtained by a fitting procedure for wurtzite 
GaN and A1N in units of Nm" 1 . 



GaN 



A1N 



ki 
k 2 

he 

k3a 



22.45 
14.96 
2.30 
0.36 



23.69 
19.55 
1.07 
-0.34 
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FIGURES 



FIG. 1. Schematic of the PFCM for wurtzite GaN or A1N along the [0001] direction. 

FIG. 2. Calculated bulk dispersion of wurtzite GaN (dotted lines) and A1N (solid lines) for 
the longitudianl phonons propagating along the [0001] direction based on the PFCM. The phonon 
symmetry of GaN at the zone-center is given. 

FIG. 3. (a) Calculated phonon dispersion of a (OOOl)-oriented (GaN)s/(AlN)5 SL for phonons 
propagating along the [0001] direction. Here, d is the period of the SL. (b) The corresponding 
atomic displacement patterns at the zone-center. The modes are numbered from the top in order 
of decreasing frequency. For clarity, the displacements are displayed perpendicular to the [0001] 
direction although they are actually pararell to the [0001] direction. The dotted line represents 
interficial N atomic plan. 

FIG. 4. Comparison of the bulk phonon dispersion of wurtzite GaN and A1N along the [0001] 
direction and the confined-mode frequencies of a (0001)-(GaN)io / (AlN)io SL mapped according to 
Eq. §). 
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Figure 1 . Zi et al. 
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